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Since the initial investigation by Matsui and Satz heavy quark bound states at finite temperature 
have been subject to numerous studies. The derivation of a finite-temperature potential from first 
principles was attempted only recently however, by generahsing the Schrodinger equation which 
is successfully employed for the description of quarkonia at zero temperature to a thermal setting. 
In this note the finite-temperature static potential is derived to leading order using resummed 
perturbation theory. The modification of the heavy quarkonium spectral function by an imagi- 
nary part of the potential appearing at finite temperature is discussed. Additionally, the extent of 
possible corrections due to non-perturbative processes is assessed by employing real-time lattice 
techniques based on kinetic theory. 
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1. Introduction 

The purpose of this brief note is to present a recently suggested definition of the finite-temperature 
static potential from first principles, which can be directly related to the quarkonium spectral func- 
tion. This so called real-time static potential |[T]] is calculated using resummed perturbation theory 
to leading order. The method of real-time lattice simulations is introduced to analyse corrections 
due to non-perturbative physics present in the classical limit of the theory. For related work on this 
subject see also [§, |3|]. 

2. Real-time static potential 

In analogy with potential models at vanishing temperature a potential V{t,r) is defined by assuming 
the time evolution of the mesonic correlator, 

iC^'{t,r)= |^/3x^,|/(^,x+^)7^WvA(f,x-|)v/(0,0)7^ 1/^(0,0)), (2.1) 

to be governed by the following Schrodinger equation, which has been shown to remain valid in a 
thermal setting at leading order [Q] : 

C2i(f,r). (2.2) 

M denotes the mass of the constituent quarks while W represents a straight Wilson line connecting 
the quark fields at time t. The notation follows the usual conventions in the Schwinger-Keldysh for- 
malism. V{t,Y) is deliberately introduced as a complex quantity with the imaginary part parametris- 
ing damping effects induced by the medium, thus generalising the concept of a potential to a ther- 
mal setting. To parametrise the evolution of the correlator at large times the potential V{y) is 
subsequently defined as the infinite time limit of y(f,r). Focusing on infinitely heavy quarks the 
amplitude C^^ (f,r) is represented by a Wilson loop of spatial extent r and temporal extent t upon 
introduction of another point splitting. Ambiguities due to time and path ordering can be avoided 
by replacing the Wightman propagator with the time-ordered propagator ' (? , r) at positive times. 
The static potential, which is referred to as the real-time static potential, is thus obtained from the 
relations 

V(r) = lim^^^^^ and /c"(?,r) ^ Ijr < #1 ^> . (2.3) 
t^oo C^^(t,Y) N \ ^ i 

The Wilson loop is assumed to be oriented in the z-direction, with ^ indicating the time-ordering. 
For convenience a temporal gauge is chosen fixing the temporal Wilson lines to unity as indicated 
by dotted lines. Summing up the series of diagrams (Figure ||), resulting from an expansion of 
the Wilson loop to first order in g^, the following expression is obtained for the real-time static 
potential, ^ 

y(r) = \imid,C''{t,r) = limlg^CF [ -^e"" {I -cos k\)iG\\{k) (2.4) 

where G^^{k) is the Fourier transform of the time-ordered HTL gluon propagator, which is related 
to the retarded propagator G^{k) and Wightman propagator G^' {k) via 

G"=G« + Gi2^GVe-'^'"G2i. (2.5) 



id,C^\t,r) 



2M-^ + V{t,r) 
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Figure 1: Symbolic representation of the Wilson loop to &{g'^). Temporal Wilson lines, corresponding to 
simple identity matrices, are represented by dotted lines. The end points of the time-ordered gluon propaga- 
tor, represented by a wiggly line, are integrated along the solid lines. 



Using the identity 

in5(co) = lim , (2.6) 

the real-time static potential is obtained from the static limits of the afore mentioned propagators, 

d^k 1— cos^^r 



y(r) 



lim «^(Gf3(«,k) + G^ (co,k)), 



(2.7) 



{27iy (/t3)2 

with a factor co^ cancelling the quadratic divergency of the temporal gauge propagator. The retarded 
propagator G'^ in temporal gauge is explicitly given by the expression 

1 kiki^ 1 



cm 



k^-ULik) w2 



The longitudinal and transversal self-energies, 



nL(^) 



(0 , G)+|k| 



1 



and Urik) 



1 



lo)\k 



k^ , (0 + \k\ 
log- 



(0 



(2.8) 



(2.9) 



develop an imaginary part for CO < \k\. The Wightman propagator is readily obtained from the 
retarded propagator by using the KMS condition G^'(^) = —2i{nB{co) + l)ImG^{k), which relates 
both propagators in thermal equilibrium {nB{(o) is the Bose distribution). An analytic expression 
for the real-time static potential is finally obtained by inserting the static limits of the propagators 



into (|2.7D, with the real part corresponding to the usual Debye-screened potential: 

d^k ,-x A 1 .%n?T. 1 



V(r) 



[In) 



[ \ — cosfc^r) 



k2 + m2 /3 |k[(k2 + m2)2 



471 



-mor 



mu + - 



dzz 



1 



sin (morz) 
niorz 



(2.10) 



27rj8 Jo (z2 + l)2 

The heavy-quarkonium spectral function p{co) is obtained by reinserting the real-time static po- 
tential into the finite-mass Schrodinger equation ( ^!2| ) and employing the relation: 



dte""'iC^^{t,0). 



(2.11) 



The spectral function for bottomonium shows the expected widening of the resonance peak with 
increasing temperature which is induced by the imaginary part of the potential (Figure^. Being the 



integral transform ( |2.7[ ) of the Wightman propagator, the imaginary part is associated with Landau 
damping. It is important to emphasise that, while the analytic continuation of the potential y(/j3,r) 
can be identified with the singlet free energy, the real-time static potential does not correspond to 
this quantity (see also [^). The static potential at finite temperature, as defined here from first 
principles, can in fact not be related in a straightforward way to an analogous quantity in euclidean 
space. 
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Figure 2: Melting of the bottomonium spectral function at finite temperature as obtained from the non- 



relativistic Schrodinger equation (2.2) using the real-time static potential. The figure is taken from 



3. Classical lattice gauge-theory simulations 

Since the introduced static potential is a genuine real-time quantity, the extent to which physical 
results are subject to non-perturbative corrections needs to be assessed employing real-time lattice 
techniques. Formally the Yang-Mills-Vlasov equations are discretised on a 3-dimensional spatial 
lattice keeping the time coordinate continouus and choosing a temporal gauge. A review of the nu- 
merical implementation of these equations is beyond the scope of this note and the discussion 
will therefore be restricted to the purely classical Yang-Mills theory. The partition function of the 
spatial lattice takes the following form, 



-I 



■Ui, 



1 



(3.1) 



£rr(i 

l'<J 

with spatial links Ui corresponding to discretised colour-magnetic fields and colour-electric fields 
defined via Ui = iEiUi. The plaquette is denoted as Uij. A discretised form of Gauss law is in- 
troduced via G{x) = Y^i \_Ei{x) — U -i{x)Ei{x — T)U^i{x)\ = 0. To evaluate the statistical expectation 
value of time dependent quantities, the partition function is supplemented by the evolution equa- 
tions 

Ui{x) = iEi{x)Ui{x) , Ei = Y^E^'T" , £f (x) = -2ImTr[T" £ Uij{x)] , (3.2) 

which follow from a variation of the discretised Yang-Mills action with respect to link variables. 
The thermalisation algorithm, generating the ensemble of configurations according to the statistical 



weight appearing in (3.1), is summarised as follows: 

1. Pre-generate the spatial gauge links Ui with a 3d Monte Carlo simulation. 



2. Generate the electric fields from a gaussian distribution [cf. eq. (3.1)]. 

3. Project onto the space of physical configurations, satisfying the Gauss law. 

4. Evolve the fields using the EOM, and repeat from step 2, until the fields have thermalised. 

For details on the implementation of the full Hard-Loop improved theory see [^. Focusing on the 
imaginary part of the potential, which remains present in the classical limit, a rectangular Wilson 
loop of spatial extent r = |r| and temporal extent t was measured using classical and HTL-improved 
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Figure 3: Time-dependence of the imaginary part of V (f , r) as obtained from lattice-regularized perturbation 
theory and classical simulations. Below is an overview of the results in the large-time hmit [3] for j3 = 16. 
Results from classical and full HTL-improved simulations {amo > 0) agree within error bars. 

lattice simulations. In Figure ^ the time dependence of the correlator is compared to a lattice- 
regularised perturbative result [^. Non-perturbative corrections were found to amplify the imag- 



inary part of the potential, which is extracted according to ( |2.3| ) in the large time limit, by up to 
100%, widening the quarkonium peak in Figure 2 but leaving the qualitative structure unchanged. 

4. Conclusions 

The non-relativistic Schrodinger equation used in potential models for quarkonia at zero tempera- 
ture is generalised to a thermal setting. The corresponding real-time static potential is derived in 
leading order perturbation theory from first principles, with an imaginary part inducing the melt- 
ing of the quarkonium spectral function expected at finite temperature. Non-perturbative processes 
present in the classical limit of finite-temperature gauge theory amplify the imaginary part of the 
potential by up to 100% and thus introduce an additional widening of the quarkonium resonance. 
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